An optimized intermolecular force field for hydrogen-bonded organic molecular crystals using atomic multipole electrostatics
نویسندگان
چکیده
We present a re-parameterization of a popular intermolecular force field for describing intermolecular interactions in the organic solid state. Specifically we optimize the performance of the exp-6 force field when used in conjunction with atomic multipole electrostatics. We also parameterize force fields that are optimized for use with multipoles derived from polarized molecular electron densities, to account for induction effects in molecular crystals. Parameterization is performed against a set of 186 experimentally determined, low-temperature crystal structures and 53 measured sublimation enthalpies of hydrogen-bonding organic molecules. The resulting force fields are tested on a validation set of 129 crystal structures and show improved reproduction of the structures and lattice energies of a range of organic molecular crystals compared with the original force field with atomic partial charge electrostatics. Unit-cell dimensions of the validation set are typically reproduced to within 3% with the re-parameterized force fields. Lattice energies, which were all included during parameterization, are systematically underestimated when compared with measured sublimation enthalpies, with mean absolute errors of between 7.4 and 9.0%.
منابع مشابه
Polarizable atomic multipole X-ray refinement: application to peptide crystals
Recent advances in computational chemistry have produced force fields based on a polarizable atomic multipole description of biomolecular electrostatics. In this work, the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field is applied to restrained refinement of molecular models against X-ray diffraction data from peptide crystals. A new formalism is also de...
متن کاملAtomic forces for geometry-dependent point multipole and Gaussian multipole models
In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise because of (1) the transfer of torque between neighboring atoms and (2) the dependence of multipole moment on...
متن کاملPolarizable Atomic Multipole Solutes in a Generalized Kirkwood Continuum.
The generalized Born (GB) model of continuum electrostatics is an analytic approximation to the Poisson equation useful for predicting the electrostatic component of the solvation free energy for solutes ranging in size from small organic molecules to large macromolecular complexes. This work presents a new continuum electrostatics model based on Kirkwood's analytic result for the electrostatic...
متن کاملSimulation of liquid imidazole using a high-rank quantum topological electrostatic potential.
Rigid body molecular dynamics simulations were carried out on pure liquid imidazole at four different temperatures and at 1 atm. Imidazole, which is important both in life science and materials science, is one of the simplest molecules to possess both a lone pair and a π system. These two features are known to benefit from multipolar electrostatics. Here the electrostatic interaction is governe...
متن کاملPolarizable Atomic Multipole-based Molecular Mechanics for Organic Molecules.
An empirical potential based on permanent atomic multipoles and atomic induced dipoles is reported for alkanes, alcohols, amines, sulfides, aldehydes, carboxylic acids, amides, aromatics and other small organic molecules. Permanent atomic multipole moments through quadrupole moments have been derived from gas phase ab initio molecular orbital calculations. The van der Waals parameters are obtai...
متن کامل